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O ! Abstract 

q ' Continuing work initiated in an earlier publication [Ichita, Yamada and Asada, Phys. Rev. D 

Q ■ 83, 084026 (2011)], we reexamine the post-Newtonian effects on Lagrange's equilateral triangular 

solution for the three-body problem. For three finite masses, it is found that a triangular config- 
q " uration satisfies the post-Newtonian equation of motion in general relativity, if and only if it has 

the relativistic corrections to each side length. This post-Newtonian configuration for three finite 
masses is not always equilateral and it recovers previous results for the restricted three-body prob- 
lem when one mass goes to zero. For the same masses and angular velocity, the post-Newtonian 
triangular configuration is always smaller than the Newtonian one. 
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I. INTRODUCTION 



One of classical problems in astronomy and physics is the three-body problem in the 
Newtonian gravity (e.g., lH3|). The gravitational three-body problem is not integrable 
by the analytical method. As particular solutions, however, Euler and Lagrange found 
a collinear solution and an equilateral triangular one, respectively. The solutions for the 
restricted three-body problem, where one of three bodies is a test mass, are known as 
Lagrange points L±, L2, £3, L4 and L5 Lagrange's equilateral triangular solution has also 
a practical importance, since it is stable for some cases. Lagrange points L 4 and L 5 for the 
Sun- Jupiter system are stable and indeed the Trojan asteroids are located there. For the 



Sun-Earth system, asteroids were also found around L 4 by the recent observations 



4]- 



Recently, Lagrange points have attracted renewed interests for relativistic astrophysics 



101, where they have d.seussed the relativistic corrections for Lagrange po,nts BB, and 
;he gravitational radiation reaction on L 4 and L 5 analytically [7fl and by numerical methods 
8MlO[. It is currently important to reexamine Lagrange points in the framework of general 
relativity. As a pioneering work, it was pointed out by Nordtvedt that the location of 
the triangular points is very sensitive to the ratio of the gravitational mass to the inertial 
one [llj. Along this course, it might be important as a gravity experiment to discuss the 
three-body coupling terms in the post-Newtonian (PN) force, because some of the terms 
are proportional to a product of three masses as Mi x M 2 x M 3 . Such a triple product 
can appear only for relativistic three (or more) body systems but cannot for a relativistic 
compact binary nor a Newtonian three-body system. 

It was shown by Ichita et al., including the present authors, that a relativistic equilateral 
triangular solution does not satisfy the equation of motion at the first post-Newtonian (1PN) 
order except for two cases 12j; (1) three finite masses are equal and (2) one mass is finite and 
the other two are zero. Hence, it is interesting to investigate what happens at the 1PN level 
for three unequal finite masses in Lagrange's equilibrium configuration. For the restricted 
three-body problem, on the other hand, Krefetz found a relativistic triangular solution by 
adding the corrections to the position of the third body j^]. For three general finite masses, 
we shall look for a relativistic equilibrium solution that corresponds to Lagrange's equilateral 
triangular one. 

Throughout this paper, we take the units of G — c — 1. 
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II. NEWTONIAN EQUILATERAL TRIANGULAR SOLUTION 



First, we consider the Newtonian gravity among three masses denoted as Mj (1 = 1,2,3) 
in a circular motion. The location of each mass is written as r /, where we choose the origin 
of the coordinates as the common center of mass, so that 



M x n + M 2 r 2 + M 3 r 3 = 0. (1) 

We start by seeing whether the Newtonian equation of motion for each body can be satisfied 
if the configuration is an equilateral triangle. Let us put 7*12 = ^23 = r%i = a, where we 
define the relative position between masses as 

ru = ri- rj, (2) 

and rjj = \rjj\ for I, J — 1, 2, 3. Then, the equation of motion for each mass becomes 

dn M 

-d¥ = (3) 

where M denotes the total mass ^Mj. Therefore, it is possible that each body moves 
around the common center of mass with the same orbital period. Figured] shows a triangular 
configuration for general masses. 
Eq. dSJ) gives 

9 M 

-n = ^, (4) 

where o;n denotes the Newtonian angular velocity. The orbital radius 77 = \rj\ of each body 
with respect to the common center of mass is obtained as 



n = aJ v\ + u 2 v 3 + J^i ( 5 ) 



r 2 = aJuf + u t u 3 + v\, (6) 



r 3 = ay v\ + v x v 2 + v\, (7) 
where we define the mass ratio as vi = Mj/M. 

III. POST-NEWTONIAN EQUILATERAL TRIANGULAR SOLUTION 

Next, let us study the dominant part of general relativistic effects on this solution. 
Namely, we take account of the term at the 1PN order by employing the Einstein-Infeld- 



Hoffman (EIH) equation of motion in the standard PN coordinate as 13Ml5| 
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where vj denotes the velocity of each mass in an inertial frame and we define 



(8) 



nu 



ru' 



(9) 



Note that Eq. ([8]) for the EIH equation expresses the acceleration of each mass, where the 
force exerted on one mass is divided by the mass. The PN force includes a product of three 
masses, whereas the acceleration by Eq. §E§ does that of two masses. 

We consider three masses in a circular motion with the angular velocity u, so that each 
rj can be a constant. In addition, the common center of mass remains unchanged for the 



equilateral triangular configuration as shown in 



121 ] . Hence, the PN location rj and orbital 



radius 77 of each body are unchanged from the Newtonian ones. As a consequence, the 



equation of motion for Mi can be written as 



12| 



9 M 

-u n = — T ri 

a 6 

where <5eihi denotes the PN terms defined as 

1 M 2 1 



'EIHl; 



(10) 



'EIHl 



16 a 3 + v 2 v 3 + 



X |{16(z/2 + V 2 U$ + ^f)[3 - (^1^2 + + Wl)] 

+ 9u 2 u 3 [2(u 2 + u 3 ) + + Av 2 v 3 + vl\}rii 
+ 3v / 3^ 2 ^3(^2 - z*$) (5 - 3z/i)n ± i|, 



(11) 



using tlx = T\jr\ and n±i = V\/r\UJ defined as the unit normal vector to r%. Eqs. ( fTOl) and 
( ITT]) seem to disagree with Eqs. (31) and (32) in 12J. However, it is not the case. This 
equality can be shown by noting that the angular velocity in the PN term, which appears 



at Eqs. (31) and (32) in 12], is equal to M/a 3 . One can obtain the equation of motion for 
M.2 and M3 by cyclic manipulations asl— > 2 — > 3 — > 1. These expressions show that the 
equilateral triangular solution is present at 1PN order only for two cases; (1) three finite 
masses are equal and (2) one mass is finite and the other two are zero 16]. 



IV. POST-NEWTONIAN TRIANGULAR SOLUTION FOR THREE FINITE 
MASSES 

For the restricted three-body problem, an inequilateral triangular solution was investi- 
gated {5), O]. Hence, for three finite masses, we study a PN triangular configuration. 
Let us denote each side length of a PN triangle as 



r u = a(l + e LJ ), 



(12) 



where eu denotes the non-dimensional correction at the 1PN order (see Fig. [I]). Here, if 
all the three corrections are equal (i.e. ei 2 = £23 = £31 = e), a PN configuration is still 
an equilateral triangle, though each side length is changed by a scale transformation as 
a — > a(l + e). Namely, one of the degrees of freedom for (e 12 , £23, £31) corresponds to a 
scale transformation, and this is unphysical. In other words, we take account of only the 
corrections which keep the size of the system but change its shape. For its simplicity, we 
adopt the arithmetic mean of three side lengths in order to characterize the size of the system 
as 



ri2 + r 2 3 + r 3i 



1+ g (£12 + £23 + £31 ^ 



(13) 



(see Appendix |A] for possible choices of fixing the unphysical degree of freedom). This 
arithmetic mean of the PN triangle is chosen to be the same as a side length of the Newtonian 
equilateral triangle as 



1 + g (£12 + £23 + £31) 



a. 



(14) 



so that the degree of freedom for a scale transformation can be fixed. Otherwise, a degree of 
freedom for a scale transformation would remain so that an ambiguity due to the similarity 
could enter our results. Thus, we obtain a constraint on (ei 2 , £23, £31) as 



£12 + £23 + £31 = 0. 



(15) 
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Hence, we look for the remaining two conditions for determining (512, £23, £31) in the follow- 
ing. 

We assume a circular motion of each body, where the angular velocity of each mass is 
denoted asw; (1 = 1, 2, 3). At the 1PN order, the equation of motion for Mi becomes 



-u'n =M 2 — + M 3 — + <W 

'21 '31 

M 3M 1 

= — s-n 



a 3 2 a 2 + U2 u 3 + v \ 
x {[z/ 2 (^i - v% - l)ei2 + ^3(^1 - f 3 - l)e 3 i]ni + v 7 ^^^ - £31)1*1.1} 
+ (5 EIH i + 0(2PN). (16) 

Note that each mass location rj may be different from the Newtonian one, because the 
origin of the coordinates is chosen as the common center of mass in the 1PN approximation. 
However, we can replace rj with the Newtonian location of Mi because of the following 
reason. The two terms including rj of Eq. ( fTBl) are expanded as 

-w?n = -tu 2 r m - c4r PN1 + 0(2PN), (17) 

M M M 

— ~ri = — T r m ~r PN i, (18) 

a A a 6 a 6 

respectively, where tni and rpNi denote the Newtonian location and the 1PN correction, 
respectively. By using Eq. (j3J), Eqs. ffTTj) and (ITSj) imply that the 1PN corrections to t*j 
cancel out in Eq. ( !T6|) . 

Furthermore, nx and n_i_i also have PN corrections. However, these corrections multiplied 
by £12 (or £31) make 2PN (or higher order) contributions in Eq. f)16p and hence they can 
be neglected. Also in <5eihij 1PN corrections to rii and n±i lead to 2PN, since they are 
multiplied by 1PN term as M 2 /a 3 . We obtain the equation of motion for M2 and M3 by 
cyclic manipulations asl— > 2 — > 3 — > 1. 

The PN equilibrium configurations can be present if and only if the following conditions 
(a) and (b) hold, (a) Each mass has to satisfy the EIH equation of motion and (b) a 
triangular configuration does not change with time. Condition (a) is equivalent to (a') the 



6 



coefficients of n±i in the equation of motion for each mass are zero: 

IM, w 

£12 -e 3 i-- — [y 2 - z/ 3 )(5 - 3ui) = 0, (19) 
o a 

S2z-e l2 -\ — (z/ 3 - z/i)(5-3z/ 2 ) =0, (20) 
o a 

1M. 

£31 - ^23 - o — (^1 - z/ s)(5 - 3z/ 3 ) = 0. (21) 
o a 

Condition (b) is restated as (b') the angular velocity for each mass is the same in order to 
keep the distance between masses unchanged: 

ul - uj 2 2 = 0, (22) 
cof - u j = Q. (23) 



Eqs. (1221) and (1231) are rewritten as 
3M 1 



9 „3 , /2 , „ „ r^^M -V2- l)ei2 + ^1 - ^3 - l)e 3 i] 

Z a Z/ 2 T ^2 ^3 T Z/3 

3M 1 r , 1N , 1% . 

~ o^"^~^ 1 2 M^2 - ^3 - 1)^23 + *Mf2 - ^1 - l)£l2 

2 a^j + z/ 3 z/i + z/f 
M 2 f 9 1 

H TH^T — ^^2^3[2(^2 + ^3) + ^ + 4z ^3 + fg] 

a 4 L 16 i/| + z/ 2 z/ 3 + z/| 

M 2 f 9 1 r , , 9l l 

+ — 1 TZ^T —^^3^i 2(1/3 + ^1 + ^3 + 4 ^i + ^1 ^ = 0, 24 

a 4 ^ 16 1/3 + V3V1 + v{ J 

3M 1 r t , i 

o~7^^ ; 2 M^l - ^2 - ljeia + ^3(^1 - ^3 - 1)^31 

2 a J I/j + Z/ 2 Z/ 3 + ^3 

3M 1 
2 a 3 z/ 2 + z/iz/ 2 + 

m 2 r 9 1 



-[z/i(z/ 3 - z/i - l)e 31 + 1/2(1/3 - f2 - 1)223] 



, 1 -, r . > . — f 2^3 [2(^2 + ^ 3 ) + ^2 + 4z/ 2^3 + ^4] 

a 4 ^ 16 z/| + ^2^3 + ^3 

+ — 1 7« _ 2 — 1 ~^W2 2 z/! + z/ 2 + z/ 2 + 4z/!Z/ 2 + z/ 2 ^ = 0, 25 

a 4 \\<6 v{ + v\v 2 + J 

respectively. It seems that (£12, £23, £31) do not always satisfy the above five conditions Eqs. 
(Tl9~]) - ( 1231) simultaneously. However, the number of independent conditions turns out to be 
two. The reason is as follows. By eliminating e 12 from Eqs. ([19]) and (|20|) . we obtain Eq. 
(|2ip . Moreover, the left-hand sides of Eqs. (j2l|) and (T25l) always vanish, if and only if Eqs. 
(fl9~|) and (|20|) are satisfied. These can be seen by direct calculations. 
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£12 = 


1 M 
24T [( " 2 


-vz) (5 


- 3^) - 


(^3 


-i*)(5 


- 3z/ 2 )] 


£23 — 


1 M 


-i*)(5 


- 3z/ 2 ) - 


fa 


-^2) (5 


- 3i*)] 


^31 = 


1 M 
24T [(Z/1 


-^)(5 


- 3z/ 3 ) - 


(>2 


-^3) (5 


-31*)] 



Thus, we obtain the expressions for (ei 2 , £23, £31) as 

(26) 
(27) 
(28) 

which recover previous results for the restricted three-body problem j^J. 

Substituting Eqs. ( 126]) and ( 128]) into Eq. ( TIB"]) , we obtain the angular velocity as 

f M \ 

U)i = o; N I 1 + ~ w pn J > (29) 

where 

wpn = ~Yq\- 29 ~ 14 :{viv 2 + z/ 2 z/ 3 + u 3 vi)]. (30) 

By using z/ x + z/ 2 + z/ 3 = 1, one can immediately show cjpn < 0, so that we find uj < wn for 
the same masses and a. In other words, for the same masses and angular velocity, the PN 
triangular configuration is always smaller than the Newtonian one. 

Table [I] shows the relativistic corrections of the distance between each body for Lagrange 
point L 4 (L 5 ) of Solar system. Here we choose Mi and M 2 as the Sun and each planet, 
respectively. In the case of the restricted three-body problem (z/3 — > 0), it is convenient to 
use Eqs. f fT9|) and f )20|) rather than Eqs. (126]) - ( 128]) . because it is natural to change not r 12 
but location of M3. Eq. ( 120]) implies that the correction of distance between each planet and 
Lagrange point L4 (L5) is approximately 5/16 of Schwarzschild radius of the Sun. Hence, 
we obtain the same values of this correction for Earth and Jupiter. The similar corrections 
are mentioned also in the previous paper [6j. The above PN effects, however, are so tiny 



that they could be neglected in the near-future measurements [Tj 

It is interesting to extend this 1PN work to higher PN orders for the gravitational wave 



physics (see |18l-l2lj for the equation of motion and compact binaries). 



V. CONCLUSION 

We reexamined the post-Newtonian effects on Lagrange's equilateral triangular solution 
for three-body problem. For three finite masses, it was found that a general triangular 



8 



configuration satisfies the post-Newtonian equation of motion in general relativity, if and 
only if it has the relativistic corrections to each side length. It was shown also that the 
post-Newtonian triangular configuration is always smaller than the Newtonian one for the 
same masses and angular velocity. Studying the correction to stability of this configuration 
is left as future work. 
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Appendix A: Choices of fixing the unphysical degree of freedom 

In stead of the arithmetic mean of the three side lengths, one might wish to use the 
geometric mean of them or the triangular area. Hence, let us mention briefly these cases. 
The geometric mean of three side lengths by Eq. (TT2]) is written up to the 1PN order as 



where 0(s 2 ) denotes the second order of (ei 2 , e 23 , e 3 i), namely at the 2PN order. This 
expression is identical with the arithmetic mean by Eq. ( {TBI . In addition, we can obtain the 
condition that these means are equal to a side length of the Newtonian equilateral triangle 

as 



Next, we consider a PN triangular area. A triangular area S is given by Heron's formula 

as 



(^12^23^31) 




(Al) 



£12 + £23 + £31 = 0. 



(A2) 




(A3) 



where 



s = 



rn + r 23 + r 3 i 
2 



(A4) 
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Hence, substitution of Eq. (Tl2|) into Eq. (1A3|) leads to 



S 



V3 



1 + g(£l2 + £23 + £3l] 



(A5) 



which corresponds to the triangular area of each side length Eq. (1A1I) . Therefore, a PN 
triangular area is equal to a Newtonian equilateral triangular one if and only if 



£12 + £23 + £31 = 0. 



(A6) 



This is identical with the condition Eq. ( 1A2I) . 

As a consequence, at 1PN level, the arithmetic mean of the three side lengths, the geo- 
metric mean of them, and the triangular area lead to the same characterizing the size of the 
system. 
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FIG. 1: PN triangular configuration. Each mass is located at one of the apexes, rj = |r/| (/ = 
1, 2, 3) denotes the orbital radius of each body. Each eij denotes the relativistic correction to each 
side length at the 1PN order. In the equilateral case, £12 = £23 = £31 = 0, namely, r\2 = = 
r 3 i = a. 
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TABLE I: The corrections for Lagrange point L4 (L5) of Solar system. Eqs. (|19p and (|20p are used 
for the evaluation. Here, we choose Mi and M2 as the Sun and each planet, respectively. Thus, 
f\2 = o,(l + £12) is the distance between the Sun and each planet. 



Planet Sun-L4 (L5) 


[m] Planet-L4 (L5) [m] 


Jupiter —0.353 


-923 


Earth -0.00111 


-923 



12 



